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Abstract 

We introduce the standard distribution width of fitness in the Bak-Sneppen 
evolution model. Through observing the fluctuation of this quantity in evo- 
lution, we find a different hierarchy of avalanche dynamics, wq avalanche. 
We also obtain the corresponding gap equation and the self-organized thresh- 
old w c . The critical exponents (r, 7 and p), which describe the behavior of 
avalanche size distribution, behavior of average avalanche size and behavior 
of relaxation to attractor, respectively, are calculated with numerical simula- 
tion. The exact master equation and 7 equation are derived. And the scaling 

relations are established among the critical exponents of this new avalanche. 
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Typeset using REVTeX 



*Email address: chlee@wipm.whcnc.ac.cn 



1 



I. INTRODUCTION 



The biological evolution theory of Darwin is well known to all. The biological evolution 
can be divided into two different types, one is micro-evolution, which focus on the dynamics 
of evolution at the level of a single population; the other is macro-evolution, which focus 
on the dynamics of evolution at the level of all species in ecosystem. In the last years, 
many scaling laws have been found in macroevolutionary data. The distribution N(m) of 
extinction sizes for families m decreases with m according to a pow-law: N(m) oc m~ a 
n. The distribution N(t) of genera lifetimes t follows N(t) oc t~ k The distribution 
of the number of genera N(g) with g species can be fitted with N(g) oc g~^ ||. The 
disappearance of species or genera is similar to radioactive decay, in which the amount of 
the "original element" decays exponentially with time, this phenomenon has been named as 
Red Queen effect [|J. The extinction events exhibit punctuated patterns f|. The punctuated 
equilibrium can interpret as a metastable state of the dynamics of that single species ||. 
If an ecosystem is at equilibrium, all species are stable. However, if one of them fluctuates 
with time, it will lead to the destabilization of its neighboring species. Thus, the avalanche 
fluctuation will sweep through the evolution. 

The existence of scaling laws in macro-evolution indicates that the dynamics of large- 
scale evolution is the result of a self-organized critical process. Based upon this view, Bak 
and Sneppen introduced a toy model for species evolution 0. It is defined and simulated 
as follows: (i) L species are arranged on a one- dimensional lattice with periodic boundary 
conditions, (ii) Random numbers /j, which are chosen from a uniform distribution between 
and 1, P(f), are independently assigned to each site (species), (iii) At each time step, 
locating the extremal site, i.e., the one with the smallest random number / m in, and mutating 
it and its nearest neighbors by assigning new random numbers in P(f). After enough 
mutations have occurred, the ecosystem reaches a statistically stationary state in which the 
distribution for the fitness of all species is statistic stationary unchanged. 

Since the Bak-Sneppen model being defined, its punctuated patterns and avalanche dy- 
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namics have been investigated widely. Observing the fluctuations of the smallest fitness, the 
model shows punctuated equilibrium behavior, i.e., it self-organizes into a critical state with 
intermittent coevolutionary avalanches of all size 0. Paczuski, Maslov and Bak studied 
the f avalanche (PMB avalanche), forward avalanche and backward avalanche and derived 
the exact scaling equations, they also established the relations among the critical exponents 
||. Recently, by monitoring the variations of the average fitness, the f avalanche (LC 
avalanche), is found, and the exact equations and scaling relations are obtained p-|IT|. As 
we know, the avalanche dynamics is a kind of macroscopic phenomenon in driven dissipative 
systems. The detailed dynamics of these systems sensitively depend on the initial configura- 
tions. But the distribution of avalanches, i.e., scaling law, does not depend on such details 
due to the universality of complexity []12 |. In the existent works, analyzing the smallest 



fitness is directly connected with the feature of individuals, but not directly represents the 
global feature; discussing the average fitness is directly connected with the global dynamics, 
but not directly manifests the difference in individuals. Can we find a new quantity, which 
can describe the global feature and the difference in individuals simultaneously? In this 
paper, we define such a quantity and observe the corresponding avalanche dynamics with 
this quantity. 

The outline of this paper is as follows: In the next section, the standard distribution 
width of fitness is defined, and its gap equation is obtained. Monitoring variations of the 
standard distribution width of fitness, the avalanche dynamics is anlyzed in details in section 
III. The exact mast equation is derived, and the scaling relations are established among the 
critical exponents. The last section is a brief summary. 

II. STANDARD DISTRIBUTION WIDTH OF FITNESS AND ITS GAP 

EQUATION 

In the Bak-Sneppen model, we denote the fitness of the i-th species as fa (i — 1, 2, • • •, L), 
L is the total number of species. The fitness represent population or living capability of 
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species, large fitness means immense population or great living capability, and vice versa. 
To describe the average population or average living capability of an ecosystem, Li and Cai 
defined the average fitness [[] |TT[ 

f av {s) = \ir,m. (i) 

i=l 

where fi(s) is the fitness of the i-th species at time s. Based upon the average fitness, we 
define the standard distribution width of fitness as 



w(s) 



i 



^E(/.w-/»w) 2 . (2) 

i=i 

From the above definition, we know that the standard distribution width of fitness simul- 
taneously manifests both global dynamics and difference in individuals. At the beginning 
of the evolution, i.e., s = 0, the fitnesses are uniformly distributed in (0,1). Thus for an 
infinite-size ecosystem [L — > oo), the initial average fitness / a „(0) equals 0.5, the initial 
standard distribution width w(0) is \^3/6. But for a finite-size ecosystem, f av (0) and w(0) 
fluctuate around the above values due to the finite size. Because of the absence of mutations 
or updates, the w(0) can not show the correlation between species, only with the evolution 
goes on, such correlation is being contained in w(s) gradually. 

As time increasing, f av (s) shows increasing tendency ||, however, w(s) shows decreasing 
tendency, see Fig. 1. After a long transience, i.e., s — > oo, the ecosystem self-organizes 
close to a critical state, and, f av (s) and w(s) approach close to different thresholds. When 
s,L — » oo, almost all species have fitnesses above the threshold f c = 0.66702 f7|, and the 
fitnesses are uniformly distributed in (/ c , 1). With definition of f av (s) and w(s), we obtain 
the relations among the thresholds 

^=^(Wc) = ^(l-/a,c). (3) 

here, w c and f avc are thresholds of w(s) and f av (s) respectively. Thus, w c = 0.096225 
and f avc = 0.83351. For the finite-size ecosystem, the self-organized critical thresholds will 
fluctuate around the above values, and this fluctuation decreases with the increase of L. 



Similar to the definition in Ref. ||||, we define the gap of w(s) as its envelop. It can be 
written as 

{iu(0) , 8 = 0, 

(4) 
mm(w(s)\s G [0, s\) , s > 0. 

The above definition means that the current gap of to(s) is the minimum of all w(s) (0 < 
s' < s). In the first row of Fig. 1, we show the gap W(s) as a function of s. With the 
increase of time, W(s) approaches close to the threshold gradually. By definition, when 
the gap falls to its next lower value, the separate instances are separated by avalanches. 
Avalanches correspond to plateaus in W(s) during which w(s) > W(s). A new avalanche 
is initiated when the gap falls, and the old avalanche terminates at the same time. Once a 
avalanche is over it will never affect the behavior of any subsequent avalanche. When the 
gap decreases, the probability for w(s) jump above W(s) increases, thus longer and longer 
avalanches happen typically. Finally, the ecosystem falls into statistically stationary state. 

Based upon the Ref. |§|9|], we can derive the gap equation of W(s). If the gap 
is changed from W(s) to W(s) + AW, the average number of avalanches occurred is 
N av = — AWL/(W(s) — w c ). Thus, when L — > oo, the average number of time steps 
required for increasing the gap from W(s) to W(s) + AW is As = (S) w , s \N av = 
— (S) W ( s j AWL/(W(s) — w c ), where (S) w ,, is the average size of avalanche of the cur- 
rent plateau of the gap. So the differential equation for W(s) is as follows 

dW(s) AW _W(s)-w c 

ds aw,As-*o As (S) w{s) L ' 

This equation describes the relaxation to attractor. If the law of (S) w ^ has been obtained, 
one can derive the law of relaxation to attractor from the above equation. 



III. AVALANCHE DYNAMICS 



We know that all self-organized critical system will exhibit power-law avalanche dynam- 
ics. This has been confirmed by a lot of models with self-organized criticality [fZ| |T~T|,|l~3| [3~9|] . 
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In this section, according to the definition of PMB avalanche || and LC avalanche ||, the 
definition of a new avalanche, Wq avalanche, is presented. Then its critical exponents is 
calculated with numerical simulation and its exact master equation is derived. Lastly, the 
scaling relations are established. 

The definition for wq avalanche is as follows. Similar to those used in Ref. we 
introduce an auxiliary parameter u> , where w(0) > wo > w c . Suppose that at time step 
s\, the current standard distribution width w(si) is less than wq- If the next standard 
distribution width w(si + 1) is larger than wo, a creation-annihilation branching process is 
initiated. The avalanche will continue to run unless w (s) becomes less than wq. This means, 
at time s, if all w(s) > Wq for (s x < s'< s — 1), the current avalanche continues. In terms 
of the above definition, the size of an avalanche is the number of time steps of subsequent 
punctuations above Wq. If the first appearance of w(s) < Wq after time s\ occurs at time 
si + S, the size of the current avalanche is S. 

Clearly, the above definition ensures the hierarchical structure, i.e., larger avalanches con- 
sist of smaller ones. With the decreasing of Wq, smaller avalanches combine into larger ones; 
if Wo is set as w c , the infinite-size avalanche appears. In the other case, with the increasing 
of Wo, larger avalanches split into smaller ones. Thus, the cutoff effects is unavoidable when 
Wo is not chosen as w c . Nevertheless, the same scaling laws can be obtained when wo is very 
close to the threshold w c . With numerical simulation, for the Bak-Sneppen model with 200 
species, chosen wq as 0.1100, we find the distribution of avalanche size follows the pow law 

P(S) oc S- T . (6) 

where, P(S) is the probability of avalanche of size S, t equals 1.63130 ± 0.06114, see Fig. 2. 
Because of the cutoff effects, the average avalanche size varies when wq changes. As same 
as those in Ref. 0,|H],|11|] , our numerical results show that the average size of Wq avalanche 



obeys 

(S) wo cx(w -w c )~^ (7) 
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In the above, w c is the threshold, the exponent 7 equals 2.57017 ± 0.06307, see Fig. 3. This 
law means the divergence of the average avalanche size satisfies pow law. We also find the 
relaxation to attractor abides 

(w - w c ) a s~ p . (8) 

here, s is the time steps, the exponent p equals 0.39049 ± 0.00780. Obviously, the value of 
p is very close to I/7. 

In the previous, with numerical simulation, we find some scaling laws of the w avalanche. 
Below, we will give some exact results. To describe the cascade process of smaller avalanches 
combining into larger ones when the parameter wq is changed, we present an exact master 
equation. Denoting the probability of w avalanche of size S as P(S,w ), when parameter 
changes from w into w + Aw , some w avalanches merge to larger w + Aw avalanches, 
at the same time, the probability of avalanche size distribution flows in and out. Because 
the termination of w avalanche is uncorrelated, the probability of a w avalanche merging 
to wo + Awo avalanche is proportional to — Awq/{wo — w c ). Thus the probability flowing 
out from P(S,w ) can be written as 

A out P = -\P(S, w )Aw /(w - w c ). (9) 

where, A is a constant. The probability flowing in P(S,w ) is 

A m P = -A w )P(S - S u w )Aw /{w - w c ). (10) 

Si=l 

When wq approaches to w c and Aw approaches to zero, we obtain the following master 
equation 

S— 1 

- (wo - ^) d P ^ Q Wo) = -\P(S, w ) + Xj2 P(Si, w )P(S - S u wo). (11) 

51=1 

As pointed out in the previous section, with the increasing of time, the standard distribution 
width shows decreasing tendency. So, when wo decreases, the first term on the right side of 
the master equation reflects the loss of S'-size avalanche, and the second term depicts the 
gain of S'-size avalanche. 
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Similar to Ref. ||10|| , defining the new variable u = — ln(wo — w c ), the master equation 
becomes into 

s-i 



dP ^ u) = -\ P (S, u) + \J2 P(Si, u)P(S - St, u). (12) 
Performing a Laplacian transformation for P(S,u), i.e., 

oo 

p(a,u) = J2 P (S,u)e- aS . (13) 
s=i 

Thus the master equation reads as 

^^ = Xp(a,u)[p(a,u)-l]. (14) 

oo 

When a = 0, p(a,u) = ^ P(S,u)e~ aS = 1, this is the normalization of P(S,u). While 

5=1 

a > 0, < p(a,u) < 1. Expanding both sides of the above equation as Taylor series 
throughout the neighborhood of the point a = 0, one can obtain 

= Hi-*(S) u + ^(s\-^(s\...} 

^(S) U -^ 2 (S\ + ^(S\...]. (15) 

oo 

where, (S k ) = Yl S k P(S,u). Comparing the coefficients of different powers of a in the 
s=i 

above equation, one gains an infinite series of exact equations. The first exact equations can 



be expressed as 



^\n(S) u = \. (16) 



Replacing u as — ln(u? — w c ), we obtain the 7 equation 

Noting the scaling law (S) oc (w — w c )~ 7 , inserting it into the above equation, one can 
easily gain 

7 = A. (18) 
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By far, we have obtained three critical exponents of the wq avalanche: r, 7 and p, which 
describe the behavior of avalanche size distribution, behavior of average avalanche size, and 
behavior of relaxation to attractor, respectively. As listed in Ref. ||, there are three other 
critical exponents: D, a and v. The avalanche dimension D is defined by the following 
scaling relation 

S oc R D . (19) 

where, S is the avalanche size (temporal duration), R is the spatial extent. The a and v are 
defined by the following scaling laws of cutoff. 

S co cx (w Q - w c y 1/a . (20) 

R co cx (w - w c )~ v . (21) 

The symbols S co and R co are the cutoff of the avalanche size and that of spatial extent of 
avalanche. 

Although there are six critical exponents, only two of them are independent. Choosing 
t and 7 as independent exponents, the other ones can be expressed as functions of them. 
Inserting the scaling law (S) oc (w — w c )~ 7 into the gap equation (Eq.(5)), and noting 
the scaling law (wq — w c ) a s~ p , one can find 

p = 7 _1 - (22) 



Integrating the equation {S) = f x c ° SP(S,w )dS, it is easy to obtain 

/Sco rSco 
SP(S,w )dS = J S 1 - T dS<xS 2 - T <x(w -w c )-^. (23) 

Comparing with (S) w oc (wq — w c )~' y , one can get 

Noting that the average number sites (n cov ) Wo covered by a u> avalanche scales near the 
critical point as [g] (n cov ) m cc R cov a R cc (w - w c ) _1 and S cc R D , so (n cov ) m cc S 1/D . 
Integrating the equation (n cov (S)) wo = J 1 co n cov (S)P(S,w )dS , one can obtain 



(S)) Wo cx S l / D - T dS cx (w - » c )-^. (25) 



Thus, the avalanche dimension D can be written as 

D = = J . (26) 

r + a-1 7(r-l) + 2-r 

Noting S oc R D and i? co a (w — w c ) _tJ , one can immediately obtain 

S C0 <xR° <x(w Q -w c y vD . (27) 
Comparing with S co oc (u> — w c ) _1 / CT , it is easy to know 

Up to now, the scaling relations among the six critical exponents are successfully estab- 
lished. The existence of two independent critical exponents means that there are two kinds 
uncorrelated critical behavior in wq avalanche. 

IV. SUMMARY 

In conclusion, by defining a new quantity, the standard distribution of fitness, a new 
avalanche, wq avalanche, is observed in the evolution of Bak-Sneppen model. We declare 
this avalanche is a different hierarchy of avalanche. Because, firstly, the observed quantity is 
different from those of Ref. JF||| and Ref. |9|-|TT|; secondly, the values of the critical exponents 



are also different from those found in Ref. [|7|,|| and Ref. fl-pT 



From the definition of the standard distribution width of fitness, we easily obtain the 
corresponding gap equation and the self-organized threshold. Apparently, with this defini- 
tion, the global feature and the difference among individuals are simultaneously described 
by the standard distribution width. So, for the evolutionary distribution with changeless 
or small changed average values, the standard distribution width is a good quantity to ob- 
serve the evolutionary dynamics. According to the definition of PMB avalanche and LC 
avalanche, the definition of the wo avalanche is presented. With numerical simulation, the 
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critical exponents r, 7 and p are obtained. Then the exact master equation is derived, from 
the master equation, the 7 equation is immediately obtained. At last, combining all scaling 
laws, the scaling relations are successfully established among the critical exponents. With 
these relations, we find there are only two uncorrelated critical behavior in the wo avalanche, 
for there exist only two independent critical exponents. 
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FIGURES 

FIG. 1. The standard distribution width and its gap. The above row corresponds to the 
variation of the gap W(s), the bottom row corresponds to the fluctuation of w(s). 



FIG. 2. The pow-law of the avalanche size distribution. 



FIG. 3. The pow-law of the divergence of the average avalanche size. 
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